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Abstract 

The Casimir energy corresponding to a massive scalar field with Dirichlet bound- 
ary conditions on a spherical bag is obtained. The field is considered, separately, 
inside and outside the bag. The renormalization procedure that is necessary to 
apply in each situation is studied in detail, in particular the differences occurring 
with respect to the case when the field occupies the whole space. The final result 
contains several constants that experience renormalization and can be determined 
only experimentally. The non-trivial finite parts that appear in the massive case 
are found exactly, providing a precise determination of the complete, renormalized 
zero-point energy for the first time. 
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1 Introduction 



Calculations of the Casimir energy in spherically symmetric situations have attracted the 
interest of physicists for well over forty years now. This is not strange, since in many 
contexts the inclusion of quantum fluctuations about semi-classical configurations turns 
out to be essential. Historically a first far reaching idea involving vacuum energies orig- 
inated with Casimir himself. He proposed that the force stabilizing a classical electron 
model arises from the zero-point energy of the electromagnetic field within and without a 
perfectly conducting spherical shell . Having found an attractive force between parallel 
plates due to the vacuum energy 0, the hope was that the same would occur for the 
spherically symmetric situation^]. Unfortunately, as Boyer || first showed, for this geom- 
etry the stress is repulsive 0, ||. Nowadays it is known that the Casimir energy depends 
strongly on the geometry of the space-time and on the boundary conditions imposed. 
This is still a very active field of research because a satisfactory understanding of the 
behaviour has not yet been found. For a number of results obtained in the last five years 
see for instance @-||. 

More recently the zero-point energy has got considerable attention in the context of the 
bag model an d chiral bag model [13|]-[T6|]. In those models the quarks and gluons 

are free inside the bag, but are absolutely confined to the interior of the bag. The sum of 
the mesonic, valence quark and vacuum quark contributions to the baryonic number have 
been found to be independent of the bag radius and of the pion field strength, being the 
vacuum quark contributions — which are analogue to the Casimir effect in QED — essential 
in the calculation of baryonic observables. The issue of regularization is certainly non- 
trivial. Under some circumstances, different regularization procedures can yield different 
results and real physical problems arise in the calculation of quark vacuum contributions 
to some barionic observables, as the energy itself. 

For massless fermions the zero-point energy was considered in [TBI. The mass- 
less fermionic field inside and outside the spherical bag was analysed in |18j. In the last 
case, a cancellation of divergences between the inner and outer spaces occurs and finite 
zero-point energies are found. Considering only the inner space, divergences appear and 
it is necessary to introduce contact terms and perform a renormalization of their cou- 
pling. Results for the massive fermionic fields contain new ultraviolet-divergent terms in 



addition to those occurring in the massless CcLS6 ; clS has been discussed in 111 911. Further 
considerations, especially on the renormalization procedure necessary in order to carry 
out these calculations, and on its precise interpretation, can be found in . 

In nearly all of the mentioned works the authors have used a Green's function approach 
in order to calculate the zero-point energy. An exception is Ref. [2(J, where, in the general 
setting of an ultrastatic spacetime with or without boundaries, a systematic procedure 
which makes use of zeta function regularization was developed. In this approach, a know- 
ledge of the zeta function of the operator associated with the field equation together 
with (eventually) some appropriate boundary conditions is needed. Recently, a detailed 

1 This could be expected because the Casimir force between plates was shown to be the same as the 
force resulting from the retarded (always attractive) van der Waals forces between the atoms of the plates. 
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description of how to obtain the zeta function for a massive scalar field inside a ball 
satisfying Dirichlet or Robin boundary conditions has been given by some of the authors 



21] , |22|| of the present work. An analytical continuation to the whole complex plane has 
been obtained and then applied to find an arbitrary number of heat-kernel coefficients. 
In the ensuing Refs. p3|, |24| the functional determinant has been considered too and, 



furthermore, the method has been also applied to spinors |25| and p-forms [fJT], [28| . For 
an alternative approach involving scalars and spinors see also the developments in [ 2j| . 
All the above considerations yield purely analytical and quite explicit formulas. In order 
to actually obtain values for the Casimir energy, however, a numerical integration had to 
be performed. This has been done in different cases, in particular for the massless scalar 
and electromagnetic fields p0| , partly reobtaining previous results. 



Here we will consider the massive scalar field inside and outside the bag, separately. 
We will discuss in detail the renormalization procedure which is necessary to apply in this 
situation and, specifically, the differences occurring with respect to the case when one 
assumes that the field occupies the whole space. The final result for the zero-point energy 
contains several constants which experience renormalization and whose physical values 
can be determined only experimentally. However, for the massive field — as is clear from 
dimensional grounds — non-trivial finite parts which depend on an adimensional variable 
involving the mass are present, and we will, in this case, be able to find here for the first 
time the complete, renormalized zero-point energy. 

The organization of the paper is as follows. In Sect. 2 we describe in detail the 
models considered and the regularization and renormalization procedure employed. In 
Sect. 3 we summarize briefly the formulas that are needed in the subsequent study of 
the zeta function of the problem at hand. We shall start with the scalar field inside the 
ball. Some additional considerations are necessary because the representation given in the 
previous articles [21 , |22|| was only applicable when the mass of the field is m < 1. Here 



we will derive formulas which are valid for arbitrary mass and very useful for numerical 
evaluations. All divergences and finite parts of the zero-point energy are calculated and its 
renormalization is performed with some care. The explicit dependence of the finite part 
in the mass is determined numerically. Afterwards, the exterior space is considered and a 
corresponding analysis is performed for this situation. Adding up both contributions, we 
see clearly how the divergences cancel among themselves as well as the influence of this 
cancellation on the compulsory renormalization process. Sect. 4 is devoted to conclusions. 
The appendices contain some hints and technical details that are used in the derivation 
of the zeta function for the non-zero mass case. 

2 Description of the model and its renormalization 

The physical system that we will consider consists of two parts: 

1. A classical system consisting of a spherical surface ('bag') of radius R. Its energy 
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reads: 

E dass = P V + aS + FR + k + ^, (2.1) 

where V = |tt-R 3 and S = inR 2 are the volume and surface, respectively. This 
energy is determined by the parameters p = pressure, a = surface tension, and F, 
k, and h which do not have special names. 

2. A quantum field <p(x) obeying the equation 

(□ + m 2 )p(x) = 0, (2.2) 

as well as boundary conditions on the surface. We choose Dirichlet boundary con- 
ditions for simplicity. The quantum field has a ground state energy 

e = Ij:4v ( 2 - 3 ) 

A (*) 

where the A(fc)'s are the one-particle energies with the quantum number k. 

For this system we shall consider three models, which will behave in a different way. These 
models consist of the classical part given by the surface and 

(i) the quantum field in the interior of the surface, 

(ii) the quantum field in the exterior of the surface, 

(iii) the quantum field in both regions together, 
respectively. 

The ground state energy is divergent and we shall regularize it by 

e, = \y.4)~^ 2s - ^> 2 - ( 2 - 4 ) 

A (fc) 

The one particle energies are determined by the eigenvalue equation 

(-A + m 2 )(p (k) (x) = \(k)(P(k)(x) (2.5) 
with Dirichlet boundary conditions on the surface 

<p w (x) =0. (2.6) 

|a;|=fl 

For the field in the interior the meaning of is obvious: (k) = (l,m,n), \^ m ,n) = 
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For the calculations we use the corresponding zeta function: 



C00 = £*«- 

(*) 



(2.7) 



In the interior region, we have 



C(«)(«) = £E(2/ + i)A 



1=0 n=0 



(l,m,n) ' 



(2.8) 



For the exterior zeta function &ext) ( s ) we m ust take into account that the radial quantum 
number is continuous. We have to subtract the Minkowski space contribution. This 
procedure is well known (see, for example, |Uf) and need not be repeated here. 

In the case of the third model we have simply to add the interior and the exterior zeta 
functions, namely 



C{total){s) = C(mi)(s) + C(ext)(s). 

Thus, in any case the regularized ground state energy is given by 



E, 



[mod) 







7;((mod)(s - r) fJ 2 



(2.9) 



(2.10) 



where mod means the model: int, ext or total. 

The divergent contributions of the ground state energy can be found most easily using 
standard heat kernel expansion 



/r(f) = xy 

(k) 



(47rt) 



3/2 



E 

5=0,1/2,1, 



B 3 t\ 



t 



H 



(2.11) 



by means of 



00 

C(s) = fu)I dt t*' 1 ^)' 



(2.12) 



where we have to take into account that in the presense of a boundary coefficients with 
half integer numbers are non-zero. As is known and can be easily found from the above 
formulas, the first five coefficients Bi (i = 0, |, 1, |, 2) can contribute to divergences of the 
energy resp. to the pole of the zeta functions. The corresponding contributions to the 
energy read: 



E, 



m 



O(div) 



64tt 2 
■2 



Bn 



327T 2 

1 

'32^ 



;Bo 



1 



- + In 

2 

1 + ln 

2 + In 



m 2 



m 



24tt 3 / 2 



B 1/2 



+ 



(2.13) 
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The heat kernel coefficients are well known (see for instance ||32|| ) and are, for the interior 
region, 



R (mt) _ 4 3 R (int) 

£> — -7r/t , a 1 i 2 



16 7T 



-27T 3 / 2 i?, B{ mt) = -7TR, Bffi = --7T 3 / 2 , 

1 3 3/2 6 2 315 # 



In the exterior region, we have: 



B 



(ext) 



B 



(ext) 



B 



(ext) 



-B 



(ext) 



1 3 
2' 2' 
i = 0,1,2,.... 



(2.14) 
(2.15) 



In order to perform the renormalization we choose as general scheme the following: 
all contributions of the heat kernel coefficients which can lead to divergences in some 
regularization have to be subtracted by means of a renormalization of the corresponding 
parameters in the classical part of the system. 

A remark is in order. The zeta functional regularization used here leaves the contri- 
butions of coefficients with half integer number finite in the limit s — > 0. This is a specific 
feature of this regularization and often much appreciated. However, in other regulariza- 
tions, as for example the proper time cutoff [J2U|J , these contributions are divergent. Thus, 
we have in each of the first two models five divergent contributions. In the third model 
we note that, in accordance with ( |2.9|) . 

rp(int) rp(ext) 
^(div) -T ^(div) ' 



rp(div) 
^^(total) 



and due to the known cancellation of divergent contributions, which is in fact due to 
( p. 15] ), only two of them remain. 

As physical system we consider the classical part and the ground state energy of the 
quantum field together and write for the complete energy 



E = E, 



(class) 



+ En 



(2.16) 



In this context the renormalization can be achieved by shifting the parameters in E( c i ass j 
by an amount which cancels the divergent contributions and removes completely the 
contribution of the corresponding heat kernel coefficients. In the first two models we 
have: 



P 



h 



F± 



h± 




a 



te- 



rn 

48^' 
m 



96' 



(2.17) 



where the upper sign corresponds to the first model and the lower sign to the second. In the 
third model there are only two contributions, which are divergent in some regularizations, 
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and the renormalization reads: 

m 3 m . 

^* + —, k^k--. (2.18) 

Within the zeta functional regularization this is a finite renormalization. 

After the subtraction of these contributions from E we denote it by E^ 6 ™^ = Eq—E'q 1 "^ 
and the complete energy becomes 

E = E {dass) + Et n) ■ (2.19) 

Within our renormalization scheme, we have defined a unique renormalized groundstate 
energy E^ en \ Often, the renormalization arbitrariness is removed by imposing some nor- 
malization conditions. In this case a natural candidate would be the requirement, that 
£,(ren) yg^g^gg f or ft _^ ^^jg j s fulfilled in our case. However, this requirement does 
not fix completely the renormalization in the first two models because a finite renormal- 
ization of h is still possible. 

With respect to this renormalization there is a qualitative difference between the first 
two and the third model. In the latter one some divergences have been cancelled when 
adding up the contributions from the interior and exterior regions. The corresponding 
terms in the classical energy do not need to be renormalized. In general, as mentioned in 



2Q1 and much earlier in connection with the renormalization of QED, those contributions 
which need a renormalization are not of quantum nature but have to be present in the 
classical part of the system and are to be determined from outside (like the electron 
mass and charge in QED) or by the dynamics of the classical part. For example, in the 
bag model one has to look for a minimum of the complete energy while varying these 
parameters. The contributions which do not need renormalization, like the one resulting 
from Bo, B\ or B2 in the third model, may be absent in the classical part and can be 
considered as pure quantum contributions. In this sense the third model contains only 
two classical parameters (a and k). 

Having in mind that only the complete energy has a physical meaning, a change of 
the normalization conditions, respectively a finite renormalization, would be equivalent 
to a change of the classical parameters and should not influence issues like that of finding 
a minimum of the complete energy by varying the parameters of the classical part. 

Special remarks are in order for the case of a massless quantum field. There only con- 
tributions from B2 remain (one performs the limit R — > in the regularized expressions, 
the BiS are proportional to m 4 ~ 2t ). In the third model these contributions are finite 
and can be considered as pure quantum corrections. They yield, together with the finite 
contributions (cf. the quantity N in the next section), the known 1/R contributions to 
the Casimir energy for a sphere and a massless field with various boundary conditions. 
However, in the models (i) and (ii) the B2 contribution is divergent and the corresponding 
1/R term in the energy must be considered as a classical contribution. Thus, in these 
cases the ground state energy for a massless field can be removed by a finite renormaliza- 
tion and the energy of the system remains formally the classical one. In that sense there 
is no Casimir effect. The same is true in the presence of a thick spherical shell (interior 
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and exterior region are separated by a finite distance with no quantum field), because 
here no cancellation between interior and exterior modes occurs. 

Similar remarks hold not only for a spherical bag but also for an arbitrarily shaped 
bag. For the infinitely thin bag a cancellation between interior and exterior modes does 
occur, while if it has a finite thickness this is not true anymore. 



3 Calculation of the ground state energy 



First we consider the interior case. As it is easily seen from equation ( |2.4j ), the task that 
remains for the evaluation of the zero-point energy is to perform a convenient analytical 
treatment of the zeta function ( |2.8| ). A precise way to obtain an analytic continuation of 
C(s) to s = —1/2 has been described in [21, 22 1 in detail, what allows us to be brief here. 
We may write the zeta function for the interior space in the form 



C(int)(s) 



N(int)(s) + Ms), 



(3.20) 



where the A^s are the contributions of the first five terms of the uniform asymptotic 
expansion of the modified Bessel functions as u — > oo and k — > oo, with vjk fixed. It 
is sufficient to subtract these five contributions in order to absorb all possible divergent 
contributions. A higher number of subtractions is possible in order to speed up the 
convergence of the remaining numerical expressions. We have called N the zeta function 
where all these asymptotic terms have been subtracted: 



sin tcs 



OO OO |- 

£(2, + 1) / *, (f ) 




m 



(3.21) 



where we have set t = 1/vl + x 2 , and rj = yi + x 2 + ln[a;/(l + yl + x 2 )]. In this formula 
the parameter s can be put equal to —1/2 under the integration and summation signs. 
The evaluation of N(l/2) is the remaining numerical task. The polynomials Di are 



D 1 (t) 
D 2 (t) 
Ds(t) 



a=0 
a=0 

E x 3 , a t 3+2a 



24 



a=0 



= — f- 

16 

25 3 
= 384 



-t 4 + -t 6 
8 16 

531 



(3.22) 



640 



221 
128* 



1105 
1152 



S 



and, in terms of their coefficients, Xi >a , the functions Ai(s) are given by 



A-x(s) 
A (s) 
Ms) 



^fTT E ^T-M 8 T ■ 2/ C g (2j + 2,-2; 1/2), 
20rr(s)f^ j! s+j 

£ ^(m J R)^T( S + j)( H (2j + 2,-1; 1/2), 



2r(s 



J- 



(3.23) 



TTT E St-M^ChC-I + • + 2j + 2s; 1/2) 

1 l s J j=o J- 



T (s + a + j + §) 



a=o r ( a + I 

It is easy to see that the above series are convergent for \mR\ < 1. Alternative representa- 
tions valid for arbitrary values of mR are derived in App. A. Using the above formulas, or 
alternatively, the formulas given in the appendix, one can perform the renomalization and 
calculate the renormalized ground state energy numerically. The result for the ground 
state energy is shown in Fig. 1 for R = 1 as a function of m. 

The dependence of E T ^ n on the radius for fixed mass is depicted in Fig. 2. It is seen 
that there is a maximimum for mR ~ 0.023. For very small values of the argument mR, 
the figure diverges to — oo, whereas for large values, the function goes to zero. In Fig. 2 
we have restricted the domain around the aforementioned maximum. 

The zero-point energy in the exterior of the spherical bag can be calculated in a much 
similar manner. Indeed, only a few changes are necessary. Subtracting the Minkowski 
space zeta function from the zeta function associated with the field outside the bag, the 
starting point here reads 



C«rt(s) = N ext (s) + £ {-lfMs) 



(3.24) 



with 



N ex t(s) 



Sill TVS 



X- 



7T 

d_ 

dx 



oo 

52(21 + 1) / dx 

1=0 „„T>I., 



mR/ v 



xv\- 

~r) 



m 



In K v {vx) — In 



—vr\ 



t2v(l + x 2 ) 1/4 , 



(3.25) 



+-D x {t) - ±D 2 (t) + ^D 3 (t) 



As is clear, one just needs to substitute the Bessel function K u for I u . The asymptotic 
contributions, as compared with those for the interior space, get an alternating sign coming 
from the asymptotics of the Bessel function K u , which exhibit this sign, as compared with 
those of the function I v (see Ref. p3|). By construction, N ext (s) is finite at s — —1/2. 
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Figure 1: Plot of the renormalized vacuum energy E r § n measured in units of the inverse 
of the radius. 



10 



E (ren) Figure 2 

m 



0.00 



-0.05 




-0.10 f 



-0.15 



-0.20 L — 

0.005 0.025 0.045 

R*m 



Figure 2: Plot of the renormalized vacuum energy E r Q en measured in units of the mass. 
The plot has been restricted to a domain around the maximum value. 
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Figure 3: Plot of the renormalized vacuum energy in units of the inverse of the radius. 



The results for Ai(s = —1/2) are given in App. A. They are the same as in the previous 
case. Again, the renormalized ground state energy can be calculated. The result is shown 
in Fig. 3. It is apparent that the slope is always negative and that the plot always gives 
positive values. It is clear that, had we plotted the same quantity in units of the mass, a 
curve with both these properties would have been obtained too. In particular, it would 
not have a maximum as the one observed for the interior case. 

In the case of the third model, i.e., for the quantum field extending to both regions 
altogether, we just have to add the two results above. As is shown in Fig. 4, there is an 
interval where the slope is positive. This would seem to leave open the possibility that a 
plot in units of the mass could exhibit a maximum. We have carefully investigated this 
possibility but the answer is negative. In other words, such alternative plot is always 
monotonically decreasing. 
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4 Conclusions 



In this paper we have developed a systematic approach to the calculation of the Casimir 
energy of a massive field in the presence of a spherical bag. The models corresponding to 
the quantum field confined to the interior and to the exterior space of the bag, respectively, 
have been discussed separately, and the differences in the renormalization of these models 
with respect to the case where the field is present in the whole space have been investigated 
in detail. 

Figures 1-4 show the quantum contribution to the renormalized ground state energy. 
This quantum ground state energy of the interior region exhibits a maximum for variable 
radius and fixed mass, as is clear from Fig. 2. Thus, we may say that if the bag is 
small enough, the quantum part of the vacuum energy induces an attractive force which 
is partly a result of our renormalization prescription by which the energy is normalized to 
zero in the infinite mass limit. This is the physically correct choice and, in order to deal 
with it, treatment of the massive case is necessary. The appearance of attractive forces 
in this kind of situation is not new, however, and has been found also when dealing with 
spinors JlTj and with the electromagnetic field on a dielectric cylinder fl34 |. 

Robin boundary conditions can be treated in complete analogy, as has been described 



in detail in Refs. ||22||-|j24||. Also the interior and exterior regions can be considered 
separately, and adding up the contributions coming from each region the same cancellation 
of divergences appear. The ground state energy of the electromagnetic field is the sum 
of the ground state energy of two scalar fields satisfying, respectively, Dirichlet boundary 
conditions and Robin boundary conditions (TE and TM modes). The Bi/ 2 heat kernel 
coefficient has opposite sign for Dirichlet and Robin boundary conditions, what leads to a 
partial cancellation of divergences between the TE and TM modes. Doing the same kind 
of calculation than the one presented here and taking the massless limit, previous results 
are reobtained j|, [5], what provides a further check of the procedure. 

Along the same lines, it would be interesting to perform the calculations for higher 
spin fields and to apply the results to realistic physical models, as the MIT bag model 
for instance. Furthermore, in complete analogy, the case of two spherical shells can be 
treated with our method. Also, the influence of dielectric media can be considered for the 
case of an electromagnetic field. 
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Figure 4: The renormalized vacuum energy represented in units of the inverse of the 
radius. 
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A Appendix: Representations for the asymptotic con- 
tributions inside the bag 

In this appendix we derive explicit representations of the Ai(s),i = — 1,...,3 (see Eq. 
( |3.23| )), which are valid for arbitrary mR. Let us start with A_i(s), which is actually the 
most difficult piece to treat. Instead of equation (|3.23|), one may use the representation 



21, 22 



v 



7T 



1=0 



oo 


(XV\ 2 2 


/ 






Vr) ~ m _ 


R/u 





Vl + x 2 - 1 



X 



(a.i; 



of which we need the analytical continuation to s = —1/2. With the substitution t 
[xvjR) 2 — m 2 , this expression results into the following one 
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z=o 



d/3 P~ 3/2 {e-* !+R2 " W ]» - e'^ 2 } 



where the Mellin integral representation for the single factors has been used. As we see, 
the /3-integral is well defined. Introducing a regularization parameter 5, A_i(s) can then 
be written as 



A_i(s) = lim A 1 _ 1 (s,S) + A 2 _ 1 (s,S) 

8^0 L 



(A.3) 
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Let us proceed with the remaining pieces. In A_i(s,5) one of the integrals can be done, 
yielding 
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(A.4) 
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For AL^s, 5), one gets 
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And adding up ( |A.4|) and (|A.5|) yields 
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. (A.5) 



(A.6) 



a form suited for the treatment of the angular momentum sum. 
To perform the summation over I, we will use 
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(A.7) 



with e — > understood. In order to get the Casimir energy we will need only the expansion 
of A_x(s) around s = —1/2. Using ( |A.7| ) one finds, after a lenghty calculation, 
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All other A$(s) can be treated in a much easier way. As a starting point for A (s) we 
choose [ED 



?7J 



—2s oo 



Using eq. ( |A.7|) , this yields immediately 
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For the remaining A^s we proceed in a different way. Let us explain the method using 
one of the contributions of A 1 (s), say: 
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This provides the immediate continuation of the sums to s = — 1/2. In fact, the infinite 
sum is convergent and in the Hurwitz zeta function the analytical continuation to s = 
—1/2 is easily performed. All pieces in Ai, i = 1,2,3 have a similar aspect and may be 
treated in the same way. Thus (we just write down the results) 
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This completes the list of expressions necessary for the analysis of the massive scalar field 
inside the bag with Dirichlet boundary conditions. 
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